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Получено решение задачи о воздействии подвижных нормальной и скручивающей нагрузок на бесконечно длинную цилиндри- 
ческую оболочку, находящуюся в упругом инерционном полупространстве. Полагается, что функции нагрузок могут быть разло- 
жены в ряд Фурье по угловой координате и интеграл Фурье по осевой координате. Движение оболочки описывается классиче- 
скими уравнениями теории тонких оболочек, а упругого полупространства ~ динамическими уравнениями теории упругости в 
потенциалах Ламе, для решения которых используется метод интегрального преобразования Фурье. Настоящая задача являет- 
ся модельной для расчёта напряжённо-деформированного состояния массива пород при неравенстве динамических нагрузок, 
передаваемых на каждый из рельсов, уложенных в тоннеле цилиндрической формы, или при вращательном движении очист- 
ных устройств в подземном трубопроводе. 


1. Задача о действии подвижной осесимметрич- 
ной скручивающей нагрузки на цилиндрическую 
оболочку в упругой среде рассматривалась в работе 
[1]. Движение периодической нагрузки вдоль ци- 
линдрической полости в упругом полупростран- 
стве изучалось в [ 2 ]. 

Аналогично [2] рассмотрим бесконечную круго- 
вую цилиндрическую полость радиуса К расположен- 
ную в упругом, однородном и изотропном полупро- 
странстве с параметрами Ламе Л, /л и плотностью р. В 
отличие от [ 2 ] полость подкреплена тонкой упругой 
оболочкой (в силу малости толщины оболочки пола- 
гаем, что окружающая среда контактирует с оболоч- 
кой вдоль её срединной поверхности) по внутренней 
поверхности которой с постоянноё скоростью с по- 
ступательно движутся апериодические нормальная и 
тангенциальная (скручивающая) нагрузки. Введём 
декартовую систему координат, ось I которой совпа- 
дает с осью полости, параллельной свободной от на- 
грузок плоской границе полупространства, а ось X 
перпендикулярна к этой границе: х<к , где И - расстоя- 
ние от оси полости до границы полупространства. 

Вследствие того, что рассматривается устано- 
вившийся процесс можно перейти к подвижной 
декартовой ( X. , У, г]=2-сі) или цилиндрической 
(г, Ѳ, ті=2—сі) системе координат. 

В подвижной системе координат для описания 
движения оболочки воспользуемся классическими 
уравнениями теории тонких оболочек ( 1 ), а для 
описания движения окружающей среды - динами- 
ческими уравнениями теории упругости ( 2 ): 
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где щ ѳ , % - перемещения точек срединной по- 
верхности оболочки; ф ц , ф ѳ , д,- составляющие реак- 
ции окружающей оболочку среды (при г=К д=а щ , 
Яо~(7 гѲ , Яг=(7 гг , где - компоненты тензора напря- 
жений в среде, у=77, Ѳ, г); ѵ 0 , р 0 - соответствен- 
но коэффициент Пуассона, модуль сдвига и плот- 
ность материала оболочки, /г 0 - её толщина; 
Р ѳ (Ѳ,гі), Р г (Ѳ,г /) - соответственно интенсивность 
скручивающей и нормальной нагрузки; 
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где и - вектор смещения упругой среды; М=с/с р , 
М=с/с, - числа Маха; с р с 5 - скорости распростра- 
нения волн расширения - сжатия и сдвига в среде. 

Так как граница полупространства свободна от 
нагрузок, то при х=И 

(7 = <7 =<7 =0. П,\ 

XX Ху ХГ} \Л) 


В случае жёсткого сцепления оболочки с окру- 
жающей средой 

И/ =г ѵ 2 =1 1 ,Ѳ,г. (4) 

Здесь и,, и д , и ц - компоненты вектора и. 

Выразив вектор смещения упругой среды через 
потенциалы Ламе [3] и=§гас1<р 1 +гоІ(^ 2 е 1 )+гоІ гоІ(^ 3 е,,), 
преобразуем ( 2 ) к виду 

Ѵ>,=М 2 ^% 7 = 1,2,3, (5) 

где М Х =М Р , М 2 = Мг= М„ 

Применив к (5) преобразование Фурье по 77 , на- 
ходим 

Ѵ 2 ф;- 777 2 Гф; = 0, 7 = 1,2,3. ( 6 ) 
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Здесь т 2 = 1 - М 2 , т 1 = т р , т 2 = т 3 = т , : , 


<р'(г,Ѳ,%) = | ср^г,Ѳ,г])е 


Представив компоненты напряжённо-дефор- 
мированного состояния среды через потенциалы 
Ламе и применив преобразование Фурье по 77 , 
можно получить выражения для трансформант пе- 
ремещений и- и напряжений а~ в декартовой 
( і=х,у,г / , ]=х,у,г /) и цилиндрической (/=/•, 0 , 77 , 
і=г,Ѳ,і)) системах координат как функции от ср’. 

Предположим, что скорость нагрузки меньше 
скорости распространения волн сдвига в окружаю- 
щей полость среде. В этом случае М 8 <\ 
(т 2 =т і =т ! >0), и решения уравнений ( 6 ) можно 
представить в виде 

Ф* = фМ+фМ, (7) 
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Здесь К п (к/) - функции Макдональда, к=т^, 
а пі - неизвестные функции и коэффициен- 
ты, подлежащие определению. 

Как показано в [4], представление потенциалов 
в форме (7) приводит к следующим выражениям 
для трансформант потенциалов в декартовой си- 
стеме координат: 
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Воспользуемся переписанными для трансфор- 
мант граничными условиями (3), с учётом ( 8 ). Вы- 
деляя коэффициенты при ё у( и приравнивая, в силу 
произвольности у, их нулю, получим систему трёх 
уравнений, из которой выражаем 8,і^0 через ко- 
эффициенты а„/. 

е ] ^,0 = Ъ^-е^^а л Ф А . (9) 

к-\ ^ «--со 


Вид определителя А и алгебраических дополне- 
ний Ад определён в [4]. Там же показано, что А - 
определитель Релея, который не обращается в 
ноль, если скорость бегущей нагрузки меньше ско- 
рости релеевской волны в упругом полупростран- 
стве. В этом случае условия существования преоб- 
разования Фурье выполняются, и для вычислений 
интегралов ( 8 ) можно воспользоваться одним из 
численных методов интегрирования, предвари- 
тельно определив коэффициенты а п] . 

Для дорелеевской скорости движущейся на- 
грузки соотношения ( 8 ) перепишутся в виде 
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у п=-со 


ё у 'Щ. ( 10 ) 


Для представления (р- (7) в цилиндрической си- 
стеме координат воспользуемся разложением 
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Подставляя в последнее выражение из (9) 
8,(^0, имеем 

<р]='Е( а *К п (к і г) + Ъ і І н (к ] г))е м , ( 11 ) 
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Подставляя (10) и ( 1 1) соответственно в выраже- 
ния для трансформант напряжённо-деформиро- 
ванного состояния (НДС) среды в декартовых и ци- 
линдрических координатах, получим новые выра- 
жения, где неизвестными будут только коэффици- 
енты а пі . Для определения последних воспользуемся 
граничными условиями (4), представив их в виде 

и ]\„ К = и Іі’ё = г І’ Ѳ ’ г ’ ( 12 ) 

00 

где и 0; (Ѳ,%) = | и ОІ (Ѳ ,г])е~^ п (Іг). 


Применяя к ( 1 ) преобразование Фурье по 77 и 
разложив функции перемещений точек срединной 
поверхности оболочки и нагрузку в ряды Фурье по 
Ѳ, для /7-го члена разложения получим: 

Ы 0щ +Ѵ 2 П ^ Ы йпѲ — 2Л / о‘зО М Он- = > 

Ѵ 2 «^0 м 0т, +^Ч„в -2 Іпи^ = С 0 {Р„ Ѳ -д пѲ ), (13) 

Зй'Лио», +2 + Ао„- =С 0 (Р т . -ц т . ), 

где е 1 = сс 0 — б 0 , е 2 — р 0 — е 0 , в 3 — — е 0 , 

Е 1 =ѵЛо м м’ Іо = А а 2 = 2| 2 +ѵ,/; 2 , 
І 0 =А %=Ѵ& + 2п\ УІ = % 2 (%о +п 2 ) 2 +2, 
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щ ш , Р щ - соответственно коэффициенты разложе- 

00 

ния и ы (Ѳ,Е) и Р‘(Ѳ,%)= | ^(Ѳ^е'^сігі в ряды 

—со 

Фурье по угловой координате 0 Ц=Ѳ; г, т=іу, Ѳ\ г). 
При V К ( (7 „ф, С/ „и ( СТ г ,ф р Г/„, (сТ 

Разрешая (13) относительно щ щ , и 0пѲ , щ ш , нахо- 
дим 



Здесь <5„ = {е х е 2 е г ) 2 ~(е&) 2 ~(е 2 % 2 ) 2 + 2^ 2 % 3 , 
8 ч і=(е 2 е 3 ) 2 -^ 2 ,<А=Д, 

5,3 = г А> 5 в і = А> 8 Ѳ 2 = (й]й 3 ) —(^2 , <5 вз = 

— /77. . О, — Ю 2 . 0. . — /7 1. , 

<5,-з = А^) 2 -<4 = 2л, <?2 = 2ѵ 0 ^о, <?3 = 

= ѵ 2 % 0 п, О, = ^ О л(4ѵ 0 -г 3 2 ѵ 2 ), 

А = 2 ^о АЧ - «Ч)> А = 2и(е, 2 -<? 0 2 ѵ 0 ѵ 2); 

А=0, А=А А=А; для индекс ,/= I соответ- 
ствует индексу г),і=2—Ѳ,і='і—г. 

Подставляя (14) в (12) и приравнивая коэффи- 
циенты рядов Фурье-Бесселя при е", получим бес- 
конечную систему линейных алгебраических ура- 
внений для определения коэффициентов а пі . 
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Рисунок. Изменения компонент НДС земной поверхности 
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После определения коэффициентов а ф исполь- 
зуя обратное преобразование Фурье, можно вычи- 
слить компоненты НДС среды в декартовой и ци- 
линдрической системах координат. При этом для 
вычисления интегралов Фурье можно использо- 
вать любой численный метод, если определитель 
полученной из граничных условий системы ура- 
внений не обращается в ноль. 

2. В качестве примера расчёты проводились для 
равномерно распределённых в интервале |ц|<0,2 м 
по внутренней поверхности бетонной оболочки 
осесимметричных нормальной и скручивающей 
нагрузок одинаковой интенсивности, движущихся 
со скоростью с=100 м/с в массиве алевролита при 
следующих значениях параметров: К = 1 м, И=2 м, 
/; 0 =0,02 м; ѵ 0 =О,2, ц 0 =12,110 3 МПа, р 0 =2,5-Ю’кг/м 3 , 
А= 1,688- ІО 3 МПа, ц=2,532-10 3 МПа, р=2,5-10 3 кг/м 3 . 
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Интенсивность нормальной нагрузки выбиралась 
таким образом, чтобы общая нагрузка по всей дли- 
не участка нагружения равнялась сосредоточенной 
нормальной кольцевой нагрузке р. 

На рисунке показан результат воздействия на- 
грузок на земную поверхность, где в поперечной 
плоскости ( 77 = 0 ) приведены кривые изменения пе- 
ремещений и=и х ц/р (м), и° у =и у ц/р (м) и нормаль- 
ных напряжений а уу =о у /р. 
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